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We present a general field-theoretic strategy to analyze three connected families of continuous 
phase transitions which occur in nonequilibrium steady-states. We focus on transitions taking place 
between an active state and one aborbing state, when there exist an infinite number of such absorbing 
states. In such transitions the order parameter is coupled to an auxiliary field. Three situations 
arise according to whether the auxiliary field is diffusive and conserved, static and conserved, or 
finally static and not conserved. 

PACS numbers: 



THE UBIQUITY OF ABSORBING-STATE 
TRANSITIONS 



This overview is devoted to a study of nonequilibrium 
phase transitions taking place between the active and the 
absorbing state of a system, as some control parameter 
is varied across a threshold value. Such transitions are 
encountered in a variety of fields ranging from chemical 
kinetics to the spreading of computer viruses From 
a theoretical standpoint absorbing state transitions form 
natural counterparts to equilibrium phase transitions. 
The transition rates used in the stochastic dynamics 
employed to model the physical phenomenon under 
consideration do not satisfy detailed balance (with 
respect to an a priori defined energy function) . In spite 
of this apparent freedom, the number of universality 
classes that the transition can fall into is incredibly 
small. Among known universality classes, that of 
directed percolation (DP) is by far the broadest. And 
indeed, in the absence of additional symmetries or 
conservation laws, as was conjectured twenty years ago 
by Grassberger an absorbing state transition will 
invariably fall into the DP universality class. The inter- 
est in absorbing state transitions was further enhanced 
as Dickman and coworkers ^ established a one-to-one 
correspondence with self-organized critical systems (see 
01 for a review on self-organized criticality). They were 
able to show that the scaling behavior observed there 
was entirely governed by an underlying nonequilibrium 
phase transition (which, as a side effect, somewhat 
tempers the mystics of SOC). 

The study of exactly which microscopic ingredients 
make an absorbing state transition not belong to the 
DP class has almost grown into a field of its own. It 
was early realized that if the microscopic dynamics 
possesses additional conservation laws the universality 
class of the transition could be different. Discrete 
conservation laws, such as the conservation of the parity 



of the number of particles 0, are known to be driving 
the transition to an independent universality class 
(the Parity Conserving or Voter class 0). A recent 
study attempts to provide a comprehensive table of all 
possible transitions involving the dynamics of a single 
field 0- Besides, the effect of a continuous symmetry 
was shown either to change the universality class of 
the transition 0,0] or to simply destroy its continuous 
nature |lOj |. The continuous symmetry present in the 
systems studied there arose from a local conservation law. 

An independent direction of research has focused on 
absorbing-state transitions in which the order parameter 
freezes into one among an infinity of absorbing states, 
but without any additional conservation law. The 
paradigmatic example of a system showing such a 
behavior is the pair contact process, initially introduced 
by Jensen and Dickman [ij, for which Mufioz and 
coworkers [T^ devised a convincing phenomenological 
picture that we shall later rely on for our analytic 
manipulations. The sensitivity of the DP class with 
respect to the coupling to an auxiliary field might 
actually provide a way out for explaining the difficulty 
of experimentally observing of the DP class > in spite 
of recent efforts (see jl^ and references therein). Among 
the few effects that may lead a transition to depart from 
the directed percolation class we list quenched disorder 
and the coupling to an auxiliary field. The present 
works focuses on the latter (which was first formalized in 
those terms by Grassberger, Chate and Rousseau 0|). 
We should also emphasize that the following exposition 
is a one-sided approach to those phenomena, relying 
solely on the field-theoretic approach, thus completely 
omitting the huge numerical efforts invested in the field. 



The existence of an infinite number of absorbing 
states (in the large-system limit) and the coupling to 
an auxiliary field are the common characteristics to the 
microscopic models that we wish to investigate here. 
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We will provide a full renormalization group picture of 
the phase transition at work in systems possessing an 
infinite number of absorbing states, with or without 
an additional conservation law. We shall rely on a 
combination of exact mappings and phenomcnological 
Langevin equations approach as a starting point for 
the calculations. We begin by introducing the three 
families of models the critical behavior of which we 
wish to understand. Then we sketch the field-theoretic 
stragegy that we will follow on the simplest example 
of an epidemic. As we turn our attention to more 
complex models we shall identify which are the new key 
ingredients that pose technical difhculties. By doing so 
we come across an unexpected connection between the 
last two families of models we want to consider. The 
conclusion section lists some open problems with the 
present approach and gives directions for possible future 
works. 



II. THREE FAMILIES OF PROCESSES 
A. Spreading of an epidemic (SIS) 

A model for the spreading of an epidemic which is 
well-knwon to epidemiolgists is the so-called Susceptible- 
Infected-Susceptible model (SIS). The population is di- 
vided into two groups, the susceptible A's and the in- 
fected i?'s, whose interactions are contamination of an A 
by a B upon encounter, 

A + B ^ B + B (1) 

and spontaneous recovery of an infected individual, 

B^ A (2) 

The motion of the individuals will be assumed to be 
diffusive, with diffusion constants Da and Db for the 
susceptible and the infected individuals, respectively. 
While the assumption of diffusive motion could itself 
be discussed (the very same reaction processes with 
underlying motion provided by a chaotic flow are used to 
model plancton population dynamics in the ocean [T^ ) 
we will take it for granted. The hope is that this simple 
assumption will apply to the majority of systems with 
such competitive rules. 

When the average density of individuals p, which is 
obviously conserved by the above rules, is high enough, 
the infection survives indefinitely. At low density, on the 
contrary, the epidemic becomes extinct exponentially 
fast. And there exists a critical density threshold pc 
separating the two steady-states, i.e. the active one, 
with ongoing spreading, from the absorbing one in which 
the extinct epidemic cannot be revived. In the absorbing 
states, A particles are freely diffusing. 



The critical properties of, say the density of infected 
individuals, in the vicinity or at the density threshold 
were studied both analytically H, 0, and numeri- 
cally [T^ ITsLligl . with a rich variety of behaviors. There 
is even still some controversy over the order of the transi- 
tion in low space dimensions (mean-field predicts a con- 
tinuous transition). 

It is not hard to coin a mean-field phenomcnological 
evolution equation for the order parameter field ■(/;, 

dti^ = DbAiP + [p/p, - - giV' (3) 

where gi is a coarse-grained contamination rate. From 
this equation we deduce that the density of infected 
individuals tp in the steady-state, undergoes a continuous 
transition from an active state in which, as p ^ pj, 
"0 ^ (p ~ Pc)^ : to an absorbing state as p < pc, with 
■0 = 0. At the mean- field level (3—1. Similarly the 
relaxation rate to the steady-state occurs over a typical 
time scale t ~ with z = 2 and there is typical 
correlation length ^ ^ |p — pd^"^, with nu — 1/2. At 
the critical point the density decays as V(i) ~ with 
i5 = 1. This set of mean- field exponents will accurately 
describe the scaling properties of the transition whenever 
the space dimension is high enough to allow diffusion 
to quickly homogeneize space fluctuations. However, in 
low space dimensions, where random walks are recurrent 
or weakly transient, correlations between long-lasting 
fluctuations will play a role and the scaling behavior 
will be modified. In a region of linear size £ density 
fluctuations are of order i'^^'^, and thus relative density 
fluctuations are of the order t^'^l'^ . From Q we see 
that the reaction relaxation time goes roughly inversely 
proportional to the density fluctuations, hence the time 
it takes for the reaction to equilibrate over the region of 
size I is of the order l'^^'^ . But this mean- field reasoning 
holds only if diffusion has acted fast enough so as to 
erase space fluctuations over the domain of size which 
takes a time of the order and for mean-field to hold 
we must have < t'^l'^ , that is d > 4. 

The features which make of the SIS a system whose 
behavior is not of DP type are the following: the order 
parameter field is coupled to an auxiliary conserved field 
standing for the local density fluctuation of particles 
(independently of their Aox B type). There exist an in- 
finite number of absorbing states, but admittedly all of 
them are equivalent. In ,2Qj Grassberger makes a distinc- 
tion between fully ergodic absorbing states (such as the 
ones resulting from the epidemic becoming extinct) and 
those which are frozen in. We now turn to an example 
of the latter. 

B. Fixed Energy Sandpiles (FES) 

The second family of systems that we would like our 
analysis to encompass has, at first sight, no relationship 
with the above. Stochastic fixed energy sandpiles (FES), 
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known as the Manna model 0, 121I |2^ . are defined as 
follows. Grains of sand are distributed on the sites of 
a lattice. Whenever a site is occupied by more than 
two grains, these are randomly distributed to the nearest 
neighbors. After a while, at sufhciently low grain den- 
sity, all the sites will have less than two grains and the 
activity remains frozen in forever. At high grain den- 
sity the toppling processes will keep occurring through- 
out the system, thus leaving a finite density of active 
-or toppling- sites in the infinite time limit. Of course, 
the system is supposed to be closed (the grains cannot 
escape though the boundaries). Identify now an active 
site with an infected individual B and a still site with 
a susceptible A. Only the S's are now diffusing, but 
otherwise the reaction rules (|1I2(I continue to roughly de- 
scribe what is taking place in the system. Interestingly, 
and this is the important contribution of Vespignani and 
coworkers, a very large number of cellular automata can 
effectively be described by the SIS model in which the 
susceptibles are static. While this slight variation on the 
model present in the previous paragraph sounds innocu- 
ous, it will prove a dramatic change for te scaling prop- 
erties. To get a first feeling of why it is so, we note that, 
as previously, the order parameter is coupled to a con- 
served field. But now this field is static, or rather its 
dynamics is slaved to that of the order parameter. We 
may expect, thus, that the fluctuations of this field will 
introduce memory effects, and lead to increased devia- 
tions to mean-field behavior. Again there exist an in- 
finite number of absorbing states, but we do not know 
a priori what their distribution is. Pastor-Satorras and 
Vespignani [2^ performed the first study of the epidemic 
with static healthy individuals, after the general charac- 
teristics of systems described by those schematic reaction 
rules had been identified by Rossi et al. by studying 
two distinct automata, the Activated Random Walkers 
and the Conserved Transfer Threshold Process. 



C. Pair Contact Process (PCP) 

The third and final family of models we will bend over 
is embodied by the Pair Contact Process (PCP), a reac- 
tion involving a single species of static particles denoted 
by C These excluding particles may either annihilate 
when on neighboring sites according to C + C ^ or 
produce an offspring C + C^C-\-C + C. But they do 
not diffuse. The order parameter of the transition is the 
density of pairs of nearest neighbor particles. At high 
branching rate the steady-state exhibits a finite density 
of particles, while at low branching the system eventually 
settles in a frozen state in which all particles are isolated. 
The order parameter of the transition is the local pair 
density. The connection between the PCP and the above 
two processes lies in the identification of a pair with a 
B particle, and of an isolated particle with an A. Un- 
fortunately there is no ready-to-use one-to-one mapping 
between a configuration of A's and i?'s and a configura- 



tion of the C's. One can draw a heuristic correspondence 
based on physical intuition, which can be written in a 
rather loose notation in a form reminiscent of the previ- 
ous paragraphs: 



A + B B + B, B 



B^A 



(4) 



While there is obviously no particle conservation one may 
easily identify a branching process similar to the contam- 
ination one and two growth-limiting processes. Note also 
that the branching of a pair induces effective diffusion for 
the i?'s. A large number of systems, such as the dimer 
reaction pj| , the dimer-dimer |25j | and dimer-trimer 
reactions, or the threshold transfer process |23| are mod- 
eled by the above mechanism. 



III. FIELD-THEORETIC TECHNIQUES: A 
UNIFIED TREATMENT 



A. Action for the epidemics 

There are well-established techniques for map- 
ping a reaction-diffusion process in which particles dif- 
fuse, branch, annihilate, and possibly exclude each other, 
onto a field-theory that exactly encodes the microscopic 
dynamics. We will not review how such a mapping is 
achieved. It suffices to know that the resulting field the- 
ory features, for each species of particles, a density field, 
the first moment of which gives the local density of par- 
ticles, and a response field (usually bearing an overbar) 
with no straightforward physical interpretation. Denot- 
ing by -0 the order parameter field and by the local 
density fluctuation we have the following action 



S[^p, V', 0] = / d'^xdt 



^{dt + (T-DBA)^ 



(5) 



where we have omitted higher order terms in the fields 
and have not specified any names for the interaction ver- 
tices. Note that the free field (f>, 4> could be integrated 
out to yield an effective action for the order parameter 
Tp and its response field ip alone. The renormalization 
group analysis of (O was performed |3, 0, 0| with the 
following results : the coupling of the order parameter to 
an auxiliary diffusive and conserved field drives the tran- 
sition to two universality classes different from the DP 
one according to whether < Da < Db or Da = Db- 
Even more surprisingly, the continuous phase transition 
is outweighed by a first order one below the upper critical 
dimension when < Db < Da ■ While a first-order tran- 
sition was indeed observed in a two-dimensional syst em, 
it has been both ruled out 0, and confirmed [s^ 
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in one dimension. In the following table we recall the 
results for the critical exponents given to first order in 
e ~ A — d {a. star superscript indicates a result holding to 
all orders in s) and provide for comparison the directed 
percolation expression ^3lJ. The exponents 77 and fj are 
the anomalous dimensions of the ip and ip fields, respec- 
tively, defined as ipR ~ ^ ^ a-nd i/iR ~ i ^, where i 
is a length scale. 



exponent 


Da = Db 


l^A < Db 


DP 


z 


2* 


2* 


^ 12 




^ 2 


^ 2 


2 - - 


V 


e 


0* 


e 
fi 


V 


e 


-0.313£ 


e 
fi 


P = v^ 


1 - ^ 
8 


1* 





(6) 



It is worth commenting on Table several exponents 
appear to be given to all orders in e. That the z exponent 
takes the value 2 is understable since the conserved and 
freely diffusive field 4> imposes its (slow) relaxation scale. 
Hence the superdiffusive behavior observed in DP (z < 2 
) when it is not coupled to any auxiliary field is overcome 
by the slowest relaxing modes, which happen to be freely 
diffusive modes (with diffusion constant Da)- This is a 
very robust property that can be seen to hold even if we 
took the omitted quartic terms in ((SJ into account. It is 
rooted in the absence of renormalization of the (j) propa- 
gator. The exactness of the correlation length exponent 
V — 2/d has a very different origin. A shift of the den- 
sity field (/) by a constant has the same effect as shifting 
the if) field mass term. This implies a exact identity be- 
tween vertex function, which carries over to the renormal- 
ized quantities, thus yielding an additional relationship 
between exponents, namely that equals the scaling 
dimension of 0, which is d/2. Unfortunately this contin- 
uous symmetry is broken by quartic terms whose effect 
in low-dimensions is ill-controlcd. This has led to some 
debate in the literature [l8|, |3^ |33| . Quite understanbly, 
the exact knowledge of some exponents (if blindly taken 
for granted) eases numerical analysis and yields higher 
precision results on the remaining exponents. 



B. One-loop-expansion for FES 

One crucial difference between the SIS and the FES is 
that the A particles become static, which can be achieved 
simply by setting Da = in ((S)). Note that an immedi- 
ate consequence of setting the auxiliary field diffusion 
constant to zero is that it ceases to impose its own dy- 
namic scale (in particular, we should now expect that the 
order parameter field V' will indeed exhibit superdiffusive 
behavior, and that will follow). The total density field 
(f) has its dynamics slaved to the local fluctuations in the 
order parameter ip, as can be see by writing the equation 
of state for 6: 



This means that the modes of describing the short-scale 
fluctuations verify 0(q, t) ~ — /* dTV'(q, t). When in- 
tegrating out the 4), (f> fields from the action Q one is left 
with effective interaction of the form 



54 / d''xd#V(x, i) / dt'AxV^x, i') 



(8) 



But at small distances (large q's) there is no difference 
between the (74 vertex and 



53 / d''a;d#V(x,t) / dtV(x,t') 



(9) 



Thus it is not surprising that a (73-like vertex is gener- 
ated already at one loop, as can be seen on the Feynman 
diagram depicted in Fig. The fact that such terms 




FIG. 1: The one-loop diagram combines the 174 and g2 vertices 
to yield effective 33 and gi couplings. A black dot on a leg 
means that the vertex is proportional to the square of the 
momentum flowing through the leg. The dotted line means 
that the danghng leg is at time earlier than the other two 
making up the vertex. 

are generated would seem to be in conflict with the local 
conservation law. But there is no contradiction as other 
terms are generated, which will turn out to be irrelevant 
from the renormalization group point of view, yet they 
will take care of preserving this essential property. We 
now rewrite the effective action for the order parameter 
field only: 



d'^xdt i){dt + At - AA)V' 

+ J dtV(x,t') 
-gii'i^ I dt'A7^(x,i') 



(10) 



gsV^V- / dtVV'(x,i') 



dt{<i>)=DB^{lp) 



(7) 



The gi and §2 vertices alone make up the directed perco- 
lation action. In fact, in all members of the FES family, 

vertices of the form Jd xdt ijjip J dt'ip(t') are gen- 
erated at one loop. To illustrate this state of affairs, 
in Fig. 121 we draw the one-loop diagram generating the 
n — 2 term. In the next section it is explained how to 
make physical sense of the action (|1U|I and to extract 
physical quantities from it. 
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FIG. 2: Two ^3 vertices are combined with a t/s vertex to 
yield an effective J d''a;di tptp ^J*^ dt'i(){t')\ interaction. 



C. The special PCP case 

In the PCP the microscopic dynamics is quite differ- 
ent. Following the suggestion of Muhoz [l^, it is not 
hard to see that the directed percolation action has to be 
supplemented by an additional interaction term of the 
form 



6Spcp [V', ^] 



d'^xdti'ipix, t) 



1 



(11) 

Note that subtracting —1 in the bracket in the right- 
hand-side of amounts to shifting the mass term of 
the field theory (which now vanishes at the mean-field 
critical point). The mass term being positive, this means 
that we are working in the absorbing phase. In other 
words, the order parameter relaxes exponentially fast to 
zero. We choose to perform our analysis with a posi- 
tive mass, and we shall carry the necessary renormaliza- 
tions in the absorbing phase, as is done in all previous 
renormalization group studies of absorbing state transi- 
tions. Working in the active phase would lead to the 
same renormalization factors, but would imply following 
a much more complicated path with exactly the same 
outcome. 

One could believe that the present exponential term 
will not contribute to determining the anomalous scal- 
ing properties of the phase transition at work by arguing 
as follows. The order parameter will tend to a constant 
value and hence this term will be suppressed exponen- 
tially fast. But this reasoning involves fixing the mass 
and letting time to infinity. We are interested in the re- 
verse limit in which time is much shorter than any typical 
critical time so that the critical fluctuations can develop 
(but of course much larger than microscopic time scales). 
And those limits do not commute. All the terms mak- 
ing up the series expansion of ^S'pcp are relevant at the 
Gaussian fixed point (the n's power of A has bare dimen- 
sion 2 -|- ne, which indicates that those terms will con- 
stitute relevant perturbation of the directed percolation 
0{e) fixed point. 

The nature of short-time and short-distance singu- 
larities (UV divergencies) dictates the universality class 
that a phase transition falls into. This is quite counter- 
intuitive since a phase transition is a large-scale collec- 
tive phenomenon. And indeed renormalization is not but 



a refined coarse-graining procedure which integrates out 
the short-time and short-distance degrees of freedom, and 
which will carry these pieces of information over to the 
macroscopic degrees of freedom (the global order param- 
eter, for instance). 



D. Analysis 

In order to study UV divergencies we follow the strat- 
egy outlines in and we expand (55pcp in powers of 
A. Then a quick inspection at the one-loop graphs tells 
us that the new upper critical dimension is dc = 6, be- 
low which the theory is super-renormalizable, and above 
which mean- field applies. From here on the notation 
£ stands for the deviation of the space dimension with 
respect to the new upper critical dimension dc = 6: 
£ = 6 — d. It is instructive to write JSpcp in the form 

dSpcAi^.i)] =A j d'^xdt^PiP{^,t) dtV(x,i') 



d'^a;d#V'(x,t)x 



g-Aj„*dt'v(x.t') / dtV(x,i') 

Jo 

(12) 



because now the brackets in the right-hand-side contain 
terms which are irrelevant in the renormalization group 
sense. A word of caution is needed here: the unfortunate 
wording irrelevant does not mean those terms are 
irrelevant physics-wise, and indeed by throwing them 
away we would simply lose the physical mechanism at 
work in the PCP (there would simply be no more phase 
transition). This merely expresses that those additional 
terms do not introduce corrections to the large scale 
effective couplings. There are many instances of systems 
described by a field theory in which RG-wise irrelevant 
terms govern the phase diagram, but not the anomalous 
scaling properties, both in equilibrium (see the review 
by Amit and Peliti |35j| ') and out of equilibrium (see 
Janssen and Schmittmann 36] for the example of a 
driven diffusive system or Deloubriere al. |37l |. more 
recently, for providing the analysis of the Pair Contact 
Process with Diffusion) . All powers of A higher than two 
therefore constitute dangerously irrelevant terms. 

However, having in mind that now the upper critical 
dimension is shifted up to dc — 6, we see that also gi 
is a dangerously irrelevant coupling to the extent that, 
as shown by Muhoz and coworkers it eventually 

controls the dynamics of the order parameter by playing 
the role of the leading nonlinear growth- limiting coupling 
(at late times). 

There are now two separate problems. The first one is 
to renormalize the field-theory with the g^ijjip /* dt''ip{t') 
vertex. But this has already being done almost twenty 
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years ago by Janssen [sl] and Grassberger and Cardy |3^ 
(this is the same field theory that describes dynamical 
percolation). And then one must compute the scaling 
dimension of the dangerously irrelevant operators at 
the dynamical percolation fixed point. The second 
issue is how to extract the scaling behavior of physical 
observables knowing that the usual scaling assumption 
breaks down. Again it may look surprising that the 
universality class is that of dynamical percolation, a 
process with no absorbing state transition. This is 
because in the POP the phase transition is driven by 
irrelevant couplings (which are absent in dynamical 
percolation). 

The scaling dimension of gi at the dynamical percola- 
tion fixed point is found to be = — 2 — |. At scale 
b ^ 1, the effective gi (b) behaves as . And then we 
use the mean-field scaling function's dependence in gi 
(which depends on 171 as I/51) which is correct to first 
order in e = 6 — d. The conclusion is that 

{^) ^b''^T{b'/''\a\,b~%g^{b)) (13) 

with ^{x, y, z) ^oc^ ^i^-v) leading order in e. We can 
therefore extract the following critical exponents 

/3 = 1-Ae, S=l-ls (14) 

which are valid to first order in e. We have used the dy- 
namical percolation expressions for the correlation time 
and correlation length exponents, z = 2 — | and = 

2 5-P 

z 21^- 

IV. CONCLUSIONS 

One of the major drawbacks of the analysis presented 
here, as far as the FES and PCP-like systems are 
concerned, is the absence of agreement with numer- 
ical simulations. That there is disagreement with 



the dynamical percolation picture in one-dimensional 
systems is not too surprising owing to the triviality 
of percolation in d = 1. This signals that there is 
very likely an intermediate dimension below which 
the universality class that we have described shifts to 
some other class. It could even be that FES and PCP 
systems behave differently in low space dimensions. 
As usual with expansions in the vicinity of the upper 
critical dimension, few results lend themselves to direct 
comparison with the simulations. We shall not try to 
review the numerical status of the systems considered 
throughout this paper, but we mention that the general 
belief is that the upper critical dimension, both for the 
FES and PCP-like systems is dc — 4, and not dc = 6 as 
advocated here. We refer the interested reader to i4(M. 



We would like to mention a list of interesting problems. 
In all three families of processes that we have consid- 
ered there is competition between a one-particle branch- 
ing process and nonlinear growth limiting processes. It 
was recently argued that binary spreading processes (first 
introduced by |4l|). that is reactions in which branching 
only occurs by pairs, would lead absorbing state transi- 
tions to a universality class different from that of directed 
percolation. While there is not even any consensus on 
that, we speculate that it is indeed so 37j], at least in the 
case of the Pair Contact Process with Diffusion (PCPD). 
In the PCPD there are two absorbing states (degenerate 
in the large size limit), but coupling it to an auxiliary 
field, possibly conserved and/or static, might reveal a 
rich variety of behaviors We leave this and other 

issues for future work. 
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